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Abstract. This paper extends the derivation of the Lagrangian averaged Euler (LAE-cr) equa- 
tions to the case of barotropic compressible flows. The aim of Lagrangian averaging is to regularize 
the compressible Euler equations by adding dispersion instead of artificial viscosity. Along the way, 
the derivation of the isotropic and anisotropic LAE-a equations is simplified and clarified. 

The derivation in this paper involves averaging over a tube of trajectories rf centered around 
a given Lagrangian flow n. With this tube framework, the Lagrangian averaged Euler (LAE-«) 
equations are derived by following a simple procedure: start with a given action, Taylor expand in 
terms of small-scale fluid fluctuations £, truncate, average, and then model those terms that are 
nonlinear functions of £. Closure of the equations is provided through the use of flow rules, which 
prescribe the evolution of the fluctuations along the mean flow. 
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1. Introduction. 

Historical Remarks. The incompressible case will be discussed first. The La- 
grangian averaged Euler (LAE-a) equations for average inc ompressible ideal fluid mo- 
tion first appeared in the context of averaged fluid models in lHolm. Marsden. and Rating 
|l998alb| . Dissipation was added later to produce the Lagrangian averaged Navier- 
Stokes (LANS-a) equations, also known as the Navier-Stokes-a equations 1 . 

Remarkably, the LAE-a equations ar e mathematically identical to the inviscid 
second grade fluid equations introduced in lRivlin and Ericksonl 195fij . except for the 
fact that the parameter a is interpreted differently in the two theories. In the case of 
LAE-a and LANS-a, the parameter a is a spatial scale below which rapid fluctuations 
are smoothed by linear and nonline ar dispersion. 

As in, for example, the work of IWhithaml |l974| on nonlinear waves, the distinc- 
tive feature of the Lagrangian averaging approach is that averaging is carried out at 
the level of the variational principle and not at the level of the Euler or Navier-Stokes 
equations, which is the traditional averaging or filtering approach used for both the 
Reynolds averaged Navier-Stokes (RANS) and the large eddy simulation (LES) mod- 
els. As such, the variational procedure does not add any artificial viscosity, a physical 
reason to consider the LAE-a or LANS-a equations as good models for incompressible 
turbulent flow. Moreover, it has been proven that the a models are computationally 



turbulent now. Moreover, it has been proven that the a n 
very attractive (see lChen et al.l |l999) : Irvlohseni et alJ |200 
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Although sharing the same general technique (use of averaging and asymptotic 
methods in the variational formulation), several alternative derivations of incom- 
pressi b le LA E-q equations exist in the literature. One of these derivations (see 
Holml Il9 99n uses the gener alized Lagrangian mean (GLM) theory developed in 



Andrews and Mclntvrel 197- 



An alter native derivation of the inco mpressible LAE-a and LANS-a equations 
was given in lMarsden and Shkollerl |2003j by using an ensemble average over the set 
of solutions of the Euler equations with initial data in a phase-space ball of radius 
a, while treating the dissipative term via stochastic variations. The derivation also 
uses a turbulence closure that is based on the Lagrangian fluct uations, na mely a 
generalization of the frozen turbulence hypothesis of Taylor (see iTavlorl |l938j ) . 

Rigorous analysis aimed at proving global well-posedness and regularity of the 
three-di mensional isotropic and aniso t ropic LANS-a equations can be fo und in, for 
example, iFoias. Holm and Titil |2002j ; iMarsden and Shkollerl |200ll l2003j . However, 
global existence for the inviscid three-dimensional Lagrangian averaged Euler (LAE- 
a) remains an open problem. 

From a computational viewpoint, numerical simulations of the a models (see 
IChen et alJ |l999j : iMohseni et al.l |2f)f)3j ) show that the LANS-a equations give com- 
parable computational savings as LES models for forced and decaying turbulent flows 
in periodic domains. For wall-bounded flows, it is expected that either the anisotropic 
model or a model with varying a needs to be used; the computational efficacy of these 
methods on such flows remains to be demonstrated. 

A s far as the compressible case is concerned, the only reference we know of is 

iHolml |2002aj . We shall discuss the relation between the work in this reference and 
the present paper below. 

We refer the interested reader to lMarsden and Shkoller fo r a more de- 

tailed history of the PDE analysis for LAE-a and LANS-a equations and to lMohseni et al.^ 
2003J for a survey and further references about the numerical aspects of these models. 

Motivation. In compressible flows there are two major problems at higher wave 
numbers, or small scales, that require special attention. These are (a) turbulence 
for high Reynolds number flows (common with incompressible flows) and (b) strong 
shocks. In both cases the challenge lies in the appropriate representation of small 
scale effects. For turbulence, the energy cascade to smaller scales can be balanced by 
viscous dissipation, resulting in the viscous regularization of the Euler equations. 

Historically, viscous dissipation has been used to regularize shock discontinuities. 
This includes adding to the Euler equation non-physical and artificial viscous terms 
and F o urier's law for h eat transfer in the shock region (see e.g.. lLiepmann and Roshkol 
|l957| : IShapirol |l953| 1. This way, the steepening effect of the nonlinear convective 
term is balanced by dissipation. We believe that Lagrangian averaging is a reasonable 
alternative way to regularize shock waves. The net effect of Lagrangian averaging is 
to add dispersion instead of dissipation to the Euler equations; that is, one adds terms 
that redistribute energy in a nonlinear fashion. In other, rather different situations, 
the tech nique of balancing a nonlin ear convective term by dispersiv e mechanisms wa s 
used by lLax and Levermorej fl9 83l for the KdV equation and by iKawaharal |l970j ; 
iKakutanTlmcriCawaliaral |l970j for plasma flows. 

The competition between nonlincarity and dispersion has of course resulted in 
remarkable discoveries, the most famous being solitons, localized waves that collide 
elastically, suffering only a shift in phase. The robustness of solitons in overcom- 
ing strong perturbations is largely due to a balance between nonlinearity and linear 
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dispersion. Note that in Lagrangian averaging, the energy redistribution mechanism 
that is introduced is nonlinear and might yield other interesting features that warrant 
further investigation. 

Another feature of the Compressible Lagrangian Averaged Navier-Stokes-a equa- 
tions (or CLANS-a equations) is that in turbulent flows with shocks, the effect of 
shocks and turbulence are simultaneously modeled by the same technique, namely 
the Lagrangian averaging method. 

Issues Addressed in This Paper. In this paper we apply the averaged Lagrangian 
methodology to derive the isotropic and anisotropic averaged models for compressible 
Euler equations. 

One goal of this paper is to present a clear derivation of the averaged equations. 
We are particularly interested in separating the two issues of averaging and modeling. 
In the derivation, a new ensemble averaging techni que is proposed and investiga ted. 
Instead of taking clouds of initial conditions, as in lMarsden and Shkollerl |20f)3l | . we 
average over a tube of trajectories rf centered around a given Lagrangian flow r\. The 
tube is constructed by specifying the Lagrangian fluctuations £ e = rf o r\~ x at t = 
and providing a flow rule which evolves them to all later times. The choice of flow 
rule is a precise modeling assumption which brings about closure of the system. 

For the incompressible case we assume that fluctuations are Lie advected by the 
mean flow (or frozen into the mean flow as divergence- free vector fields), and we 
obtain both the isotropic and the anisotropic versions of the LAE-a equations. The 
advection hypothesis is the natural exte nsion to vecto r fields of the classical Frozen 
Turbulence Hypothesis of G.I. Taylor fsee lTavlorl |l938|) stated for scalar fluctuations. 

The second goal of this work is to extend the derivation to barot ropic compres sible 
flows. This problem has already been considered by Holm fsee iHolinl |2002aj ) in 
the context of generalized Lagrangian mean (GLM) motion. In this work, an alpha 
model appears as a G LM flu i d theor y with an appropriate Taylor hypothesis closure. 
However, even though iHolml |2002a] enumerates several frozen-in closure hypotheses, 
the averaged equations are derived only for the case when the fluctuations are parallel 
transported by the mean flow. In our work we will consider a more general advection 
hypothesis to study the compressible anisotropic case. In addition, a physically based 
new flow rule is introduced to deal with the isotropic case. 

The averaging technique consists of expanding the original Lagrangian with re- 
spect to a perturbation parameter e, truncating the expansion to 0(e 2 ) terms, and 
then taking the average. It turns out that the averaged compressible Lagrangian de- 
pends on the Lagrangian fluctuations £' o nly through th ree tensor quantities which 
are quadratic in In the terminology of IHolml |2002aj these tensors represent the 
second-order statistics of the Lagrangian fluctuations. Evolution equations for these 
tensors arc derived from a core modeling assumption: a prescribed flow rule for the 
time-evolution of the fluctuations The flow rule gives us closure, allowing us to 
apply Hamilton's principle to the averaged Lagrangian and thereby derive an equation 
for the mean velocity u. 

The organization of the rest of the paper is as follows. In 321 we describe a 
general procedure for Lagrangian ensemble averaging. This procedure is then applied 
to the action for incompressible fluids in J^to demonstrate our derivation technique. 
The general procedure is applied again in 21 this time to the more complex case of 
barotropic compressible fluids. fJSlis devoted to modeling issues; here the strategy of 
modeling the evolution of Lagrangian fluctuations £' using flow rules is discussed in 
detail. In ^ we derive the averaged equations for incompressible and compressible 
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models in both isotropic and anisotropic versions. The Appendix provides technical 
details about the fluctuation calculus used throughout the paper. 

Main Results. The main result of this paper is the derivation of compressible 
Lagrangian-averaged Euler equations with 

• anisotropic modeling of fluid fluctuations — see equations (|6.2|) . 

• isotropic modeling of fluid fluctuations — see equations (|6.4I) . 

In addition, we provide an improved derivation of the incompressible isotropic and 
anisotropic LAE-a equations. 

2. General Lagrangian Averaging. A mathematical setting for a certain class 
of compressible fluid flow problems will first be given. After describing the general 
procedure for Lagrangian averaging, the specific case of the Euler action for fluids will 
be considered. 

Let M be an open subset of M. N , representing the containing space of a fluid. 
Suppose we are given a Lagrangian for a compressible fluid, no), where ip G 

Diff(M) the space of diffeomorphisms of M, O, ip) G TDiff(M) and fi a G A N (M), the 
space of A-forms on M. Fix a time interval [0, T] and let C(Diff(M)) be the path space 
of smooth maps from [0, T] into Diff(M). Then the action S : C(Diff(M)) x A N {M) -> 
K is 

S(VifMi)= LCn(t), »)(*)> Mo) dt 
Jo 

We seek an averaged action S a (rj,no), where a is a length scale characterizing the 
coarseness of the average. Taking rj and /j,q as given, we shall describe how to compute 

Remark. It is important to emphasize that for both S and S a , r\ is merely a test 
curve. It is not an extremal of the action S. We are trying to average the action S 
itself, not any fluid dynamical PDE or the solutions of such a PDE. Our final product 
S a should not depend at all on an initial choice of the test curve r\. 

Tube Initialization. The first step is to take £ e (x,i) to be a family of diffeomor- 
phisms about the identity. That is, 

for each e > 0, f (•, t) G Diff(M) for all t, and 
at e = 0, £ £ (x, t) — x for all x, t. 

Define the vector fields £' and £" via 

d 2 



? 8e 



C and r - , , 



Use £ e to construct a tube of material deformation maps that are close to ij by letting 
ij e (X,t) = ^ e (n(X,t),t), or, written more compactly, 

V'=?°V- (2-1) 

Here, A is a material point in the reference configuration. Define the spatial veloc- 
ity by u e (x,t) = rf{(jf) {x,t),t), where rf is a given material deformation map. 
Compactly written, this reads 

u* = rfo{rf)-\ (2.2) 

The map u e is a time-dependent vector field on M, i.e. for each e > 0, and for all t, 
u e {-,t) G X(M). 
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Averaging. The existence of an averaging operation ( ■ ) will now be postulated. 
The properties this operation is required to satisfy and an example of such an opera- 
tion will be given shortly. 

Relationship Between u e andu. It is desirable to have the fluctuations £ e centered, 
on average, about the identity: (£ e (a;, t)) = x for all positions x at all times t. What 
is actually needed is that for n > 1, 



de r ' 



(2.3) 



£ = 0/ 



In other words, the n-th order fluid fluctuation vector fields should all have mean 
zero. Restricting the map to be centered about the identity means simply that the 
average will not be skewed in an arbitrary direction. From l|2.2|) and (|2.3|) one can 
derive 



(u t o^(x,t)} = u{x,t). 



(2.4) 



Equation l|2.4[) shows in which sense the average of u e is u in a Lagrangian mean theory 
defined by (rf(-,t)) = rj(-,t). T his equation is closely connect ed with the generalized 
Lagrangian-mean description of lAndrews and McTntvrel[T97^ . where the Lagrangian 
mean velocity u L and the fluctuating Eulerian velocity ifi are related in a similar way. 

Density. For the non-averaged Lagrangi an L, /in is a parameter in the sens e 
of Lagrangian semidirect p r oduct theory; see lMarsden. Ratiu. and Weinsteinl |l984j ; 
iHolm. Marsden. and Ratiul 1998b]. The physical interpretation of fio is as follows. 
Since p,Q is an AT- form on M, it can be written as 



/la = po dx l A • • • A dx N 

where po is a smooth function on M. Now po(X) is the density of the fluid at the 
material point X in the reference configuration. This is in contrast to the spatial 
density p e (x,t), which gives us the density of the fluid at the spatial point x at time 
t. Defining 



one has the relationship 



// = p e dx 1 A ■ ■ • A dx 



(?7 £ )*^o = M e 



(2.5) 



(2.6) 



Fluctuation Calculus. Because u e and p e will be expanded, the e-derivatives of u e 
and p e need to be calculated. First, define 



d_ 



and 



=o 



de 2 



(2.7) 



c=0 



By differentiating l|2.2|) . one finds expressions for u' and u" in terms of it, and 
The calculations c an be performed intrinsically using Lie derivative formulae — the 
results, as found in lMarsden and Shkollerl |2003| . are 



u" = a t e + [u, a - ^vu' ■ c - vvu(e, o- 



(2.8a) 
(2.8b) 
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In these formulas, the bracket [a;, y] = £ T y is the standard Jacobi-Lie bracket of vector 
fields on M (see, for example. lAbraham. Marsden. and Ratiul ^988]). Next, define 



d 2 

p e , and p" = — T 



P e - (2-9) 



e=0 



One obtains expressions for p' and p" in terms of p, and £" by differentiating (|2.6|) 
(see the appendix for the detailed calculations). The results are: 

P ' = -dW(pO, (2-lOa) 
p" = div(div(p£' ® £)) - div(p£")- (2.10b) 

Averaging Operation. In the above development, an averaging operation has been 
implicitly used. The properties it is required to satisfy will now be spelled out. Let 
T{Y) mean the space of smooth, real-valued functions on a manifold Y. If Y is 
infinite dimensional, then smoothness is understood in the sense of infinite dimensional 
calculus with respect to, for example, suitable Sobolev topologies. These infinite 
dimensional technicalities will not be required in any detail in this paper, and so may 
be treated formally. 

As before, the set M is the containing space of the fluid and a is a small positive 
number. Let X be an appropriately chosen space of fields, designed to model "fluid 
fluctuations," on M, and consider the space Y = [0, a] x X. Assume that there is an 
averaging operation 

{■):T{Y)^T{M) 

satisfying the following properties for f,g G ^(Y), a,b e R, ip e JF([0, a]), and 
h e F{X), 

Linearity: (af + bg) = a(f) + b(g), (2.11) 
Independence: (iph) = - (J ip{e)de^J(h), (2.12) 

Commutativity: (^J f dx^j = J (/) dx, (2-13) 

(df) = d(f), where d = d t or d = d x i . (2.14) 

Here, tph G 3-(Y) is defined as the pointwise product. Note that if ip is a constant, 
then the first and second requirements are compatible. 

For compressible flow, the space of fluid fluctuations is X = X(Af). For incom- 
pressible flow, the space of divergence-free vector fields is used instead, i.e. X = 
Xdiv(-AZ). In general, X = IidcntityA', where X is the space to which the tube maps 
£ e belong. 

Example. Let p be a probability measure on the unit sphere S in X(M), and 
define the average of a (vector-valued) function f(e,w) on [0, a] x S by 



1 f c 
a Jo 



/(e, w) dp(w) de. 



s 



One checks formally that this is an example of an averaging operation that satisfies 
the desired properties. 
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3. Incompressible Flow Revisited. Before applying the averaging technique 
to the case of compressible flow, we shall first derive averaged equations for incom- 
pressible flow, equations which have already been derived in the literature. The 
presentation given here has the advantage of being easily generalized to compressible 
flows. This advantage stems from the careful use and interpretation of modeling as- 
sumptions on the fluctuations £' — only intuitive assumptions are required regarding 
the mean behavior of the fluctuations as well as a first-order Taylor hypothesis. Fur- 
thermore, great care has been taken to separate the algebraic issues involved with the 
averaging procedure from the modeling issues. 

In the incompressible case, fluid fluctuations are modeled using the volume- 
preserving diffeomorphism group on M which is denoted by Diff vo i(M). Therefore, 
the tube construction from the previous section now reads: let £ £ (a;,i) be a family of 
volume-preserving diffeomorphisms about the identity. That is, 

for each e > 0, f(',«)e Diff vo i(M) for all t, and 
at e = 0, £ e (x, t) = x for all x, t. 

This forces £'(-, t) to be a divergence-free vector field for all t. 

Averaged Lagrangian for Incompressible Fluids. Let us start with the standard 
Lagrangian 

l(u e )= I \\\u<\\ 2 dx, (3.1) 
and expand u e in a Taylor series about it: 

u £ = u + eu' + ^e 2 u" + 0(e 3 ). (3.2) 
Substituting this expansion into H3.1(l gives 




Let l(u e ) be the truncation of I to terms of order less than e 3 . Using formulas (|2.8J) . u 1 
and u" can be rewritten in terms of u, and We do this in order to write I as a 
function only of u, and Making the substitutions and rewriting in coordinates, 

*V) = J + e(u*(d t e) + uVe^ - //TV..) + y (W i )(ftn 

- 2£'> fc ^'V + 2u',. J C h C J «' + - u^e'^' V) 

+ j ((^"V + rff'ju' - e" i «> < ) dx, (3.4) 

where the notation u*. means du l /dx J . Throughout this paper, there is an implied 
sum over repeated indices. The averaged Lagrangian for incompressible flow is now 
simply l& = (1). 
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Zero-Mean Fluctuations. Before undertaking this computation, recall from ^ 
that the fluctuation diffeomorphism maps £ £ are required to have as their average the 
identity map. This statistical assumption regarding the behavior of the fluctuations 
is the first modeling assumption: 

(O = and (£") = 0. (3.5) 

This point would not be worth belaboring except that, when combined with the 
properties of our averaging operation (|2.11I2.14[) . assumption 13.5fl forces all linear 
functions of and their derivatives to also have zero mean. Applying this fact 

to (|3.4I) causes the entire 0(e) group and the second 0(e 2 ) group (i.e. the last line of 
1)3.4(0 to vanish inside the average. 

We continue analyzing 13.4(1 : the only remaining terms are (l/2)u l u l and the first 
0(e 2 ) group. Within this 0(e 2 ) group, we integrate certain terms by parts and notice 
that all terms involving time-derivatives of £' group together: 



where D/Dt is the material derivative: 

{dt + u ■ V) . 



)- 






Dt 



D 
Dt 



(3.6) 



(3.7) 



We then simplify the remaining non-time-derivative terms from (13. 4|) . integrating by 
parts to remove second-order spatial derivatives. The final expression for the averaged 
incompressible Lagrangian is 




D^' 



Dt 



<tr(V£' ■ V£')> IN 



dx. 



(3.8) 



Modeling of £' . Immediate application of Hamilton's principle to l|3.8|l does not 
yield a closed system of equations. Namely, we have initial (t — 0) data for £' but no 
way to compute this vector field for t > 0. Our approach in what follows will be to 
write down, based on physical considerations, an evolution law, or flow rule, for 

A flow rule consists of a prescribed choice of </> in the following evolution equation 
for f: 



— = 4>{u,p,C) 



(3.9) 



Given a choice of £' at t — 0, this equation will uniquely determine £' for t > 0. Let 
us assume we have a linear flow rule, 



DC 
Dt 



(3.10) 



where il 1 - 7 is allowed to depend on u and p but not on £ £ or its derivatives. The caveat 
here is that our choice of f2 must be compatible with incompressibility; in particular, 
div£' = at t = 0, and must be chosen such that £' remains divergence free as it 
evolves. At this stage, one might raise the issue of the tube £ e and request a concrete 
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description of the whole object. Such a description is unnecessary; in order to close the 
system of evolution equations resulting from l|3.8fl , we need only describe the evolution 
of the first-order fluctuation field Now defining the Lagrangian covariance tensor 

F=(£'®0 (3-11) 
and using the linear flow rule (|3.10|) , the Lagrangian i|3.8[l can be rewritten as 



l? n (u) 







a 2 r 















2 ' 3 



dx. (3.12) 



Here we have used the fact that £' must be divergence-free. 

Advection Flow Rule. The first flow rule we shall consider results from setting 
<>" = a[y. 

<'/". (3.13) 

Using the definition of the material derivative, it is trivial to see that this flow rule is 
equivalent to Lie advection of dt£,' = —£ u £,'- This advection hypothesis is the vector 
fie ld analogue of the classical Frozen Turbulence Hypothesis of G.I. Taylor introduced 
in iTaylorl Tl9 38l. This hypothesis is widely used in the turbulence community (see 
1969j for instance for usage of this hypothesis even in the sense of Lie advection 



Cocke 



of vector fields). More recently, this generalized version of Taylor hypothesis has been 
use d to achieve turbulence closure in th e derivation of incompressible LAE-q equa tions 
Csee lMarsden and Shkolierl 120011120031 ) or in the work of Holm fsee lHolml |2002aj l on 
averaged compressible models using the generalized Lagrangian mean (GLM) theory. 

The advection flow rule l|3.13|l is perhaps the most obvious choice for that is 
compatible with incompressibility. Note that if div£' = at t = Q, then differentiating 
(|3.13|) with respect to x l yields 

auiiv^ c< °- 

Therefore, div£' = for all t > 0. Using this flow rule, both anisotropic and isotropic 
models shall be developed. For incompressible flow, no other flow rules will be con- 
sidered. 

Incompressible, Anisotropic, Inhomogeneous Flow. In this case, the flow rule is 
used to derive an evolution equation for the covariance tensor F. Time-differentiating 
pij = ^ l ^'3^ an j us i n g |j3.i3[i yields the Lie advection equation dtF = —£ U F. 
Equipped with an evolution equation for F, we can apply Hamilton's principle to 
(|3.12|) and derive a closed system with unknowns u, the average velocity, and F, the 
covariance tensor. 

Carrying this out, one finds that the anisotropic LAE-a equations are given by 
the following coupled system of equations for u and F: 

d t {l - a 2 C)u + (u ■ V)(l - a 2 C)u = - gradp, (3.14a) 

divu = 0, (3.14b) 
d t F + VF • u - F ■ Vu - V/ • F = 0, (3.14c) 
where p is the fluid pressure, and the operator C is defined by 

Cu = div[Vu • F]. (3.15) 
When a = 0, the system (|3.14ai3.14b|) reduces to the incompressible Euler equation. 
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Note. Start with the generic incompressible averaged Lagrangian (|3.12[1 and sub- 
stitute the advection flow rule (|3.13|) . Now integrate the last term by parts and use 
div£' = 0. The result is 

= f [I IMI 2 - y u ■ t VVu : F l} dx > ( 3 - 16 ) 

which is exactly the Lagrangian used in iMarsden and Shkollerl 1200.1 to derive the 
anisotropic LAE-a equations. However, in IMarsden and Shk ollerl |2003| the second- 
order Taylor hypothesis 

where the orthogonality is taken in L 2 , is necessary to achieve closure. Our choice of 
modeling assumptions rendered unnecessary any such hypothesis on the second-order 
fluctuations £ ''. Second-order Taylor hyp otheses, unlike the first-order hypothesis 
retained from IMarsden and Shkollerl |2003j . do not have much precedent in the tur- 
bulence literature, as discussed above. 

Incompressible, Isotropic, Homogeneous Fluids. To model the motion of an ap- 
proximately isotropic fluid, we take the covariance tensor F to be the identity matrix, 
i.e. 

F^ = = S ij . (3.17) 

The choice of F u = 5 lJ is a modeling assumption, and will thus only be valid for flows 
which almost preserve this property. Note that (|3.17|) is strictly inconsistent with the 
advection flow rule, and thus can only be regarded as an approximation. 

For the case of incompressible isotropic mean flow, we assume that (|3.17|1 holds; 
then differentiating this equation with respect to x k and x^ and using the fact that 
£' is divergence-free, we have 

Hence 

(tr(Vr ■ VO> = (C<.) = - (Cd") = 0, 
and the Lagrangian l|3.8|) simplifies to 

e»,,„M = / M {i|MI 2 -4( 

We emphasize that this is only an approximation, so that 

J§Uo(«) « *£(«) 

along fluid trajectories u(t) for which the covariance tensor is approximately the iden- 
tity. Now using the flow rule given by (|3.13|) , the averaged Lagrangian lf n from l|3.18|l 
becomes 

\=u i tj u%(e j e k )=ni j u% (3.19) 



Dt 



dx. 



(3.18) 
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where we have used the isotropy assumption (|3.17|) . Hence, (|3.18|) becomes 



i?(u) 



A I 



- 2 \\u\ 



yllV.II 2 



dx. 



(3.20) 



This expressi on for the averaged Lagran gian in the isotropic case is identical to the 
one derived in lMarsden and Shkollerl |200l| . Now applying either Hamilton's principle 
or Euler-Poincare theory, we obtain the standard isotropic LAE-a equations: 

d t {l - a 2 A)u + (u ■ V)(l - a 2 A)u - a 2 (Vu) T ■ Au = - gradp, (3.21a) 

divu = 0, (3.21b) 

where p is the usual fluid pressure. 

4. Averaged Lagrangian for Compressible Flow. Having understood the 
incompressible case, we now turn to the compressible case. The procedure is identical 
in all aspects except we must now keep track of density fluctuations. Start with the 
reduced Lagrangian for compressible flow: 



M 



\u'\\ 2 -W(jf) )p<dx. 



(4.1) 



The fluid is assumed to be barotropic, meaning that W, the potential energy, is a 
function only of p, the fluid density. Now expand the velocity and density in Taylor 
series 



u e = u + eu' + ^e 2 u" + 0{e 3 ) 



P = P 

and also expand the potential energy W: 



e P ' + ^p" + 0(e% 



(4.2) 



1 



W(f) = W(p) + eW'(p)p' + -^(W'ip)^ + WW) + 0(e 3 ). 
Substituting these expansions into the reduced Lagrangian gives 



J M (\M'-w<r)) P 



(u-u> - W '(p)f/)p+ ^-\\u\\ 2 -W(p)j p> 

~ ((|K|| 2 + u" ■ u) - {w"{p)p l2 + WW)) p 



(4.3) 



1 (1 



+ («•«'- W'(p)p')p' + - -\\u\\ 2 W(p) p" 



2 \2 



C(e 3 ) dx. 



This expansion is now truncated, leaving out all terms of order e 3 and higher. Denote 
the truncated Lagrangian by l(u e , p e ), and define the averaged Lagrangian l a by 



l a {u,p) = (l{u\f)). 



(4.4) 



We now outline the procedure by which we arrive at a final written expression for the 
averaged Lagrangian l a . The algebra is straightforward but tedious, so details will be 
omitted. 
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1. Use equations (|2.8|) and (|2.1U|) to rewrite l|4.3[) in terms of only u, p, and the 
fluctuations 

2. Remove two kinds of terms that vanish inside the average: 

(a) linear functions of £' or 

(b) linear functions of derivatives (either spatial or temporal) of £' or 
Note: see "Zero-Mean Fluctuations" in SJ31 for justification. 

3. Carry out the averaging operation. As in the incompressible case, the only 
quantities left inside the average should be nonlinear functions of 

The end result for the averaged Lagrangian for compressible flow is 



Cnp( W ' p) 




pW{p) + a 2 



DC 



Dt 



±w'(p) (div(pe') 2 ) - \™(p) <divdiv(p£' ® O) 



dx. 



(4.5) 



We have introduced w, the enthalpy 2 , defined by 

w(p) = W(p)+pW'(p). 



(4.6) 



5. Flow Rule Modeling. In deriving the expressions (|4.5|l and l|3.8|l for the 

averaged Lagrangians, no assumptions were made regarding how the Lagrangian fluc- 
tuations £' evolve. In this section we describe one possible strategy for modeling 
Note that such a strategy is necessary to achieve closure for the evolution equations 
associated with the Lagrangians l|4.5|l or l|3.8fl . 

Preliminary Observation. Assuming £' evolves via a linear flow rule, as in (|3.10(l . 
the vector field £' appears in the averaged Lagrangian l|4.5|l only as part of the following 
three expressions 3 : 



re 



11 = {CC 



(5.1a) 
(5.1b) 
(5.1c) 



Note that F is the same Lagrangian covariance tensor from the incompressible deriva- 
tion. In terms of these quantities, the averaged compressible Lagrangian is given in 
coordinates by 



comp 



M 



1 i i 

-pu l u l 



pW(p) 



2 F 



(5.2) 



dx. 



2 Any function w satisfying Vui = (Vp)/p, where p is pressure, is called enthalpy. Our definition 
of to implies w t i = 2W'(p)p,i + pW"(p)pj = (p 2 W'(p)),i/p = p,i/p as required. 



J Similar tensors appear in H olml [1 
Lagrangian fluctuations 



; they are referred to as second-order statistics of the 
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Time-differentiating (|5.1al5.1c|l and using the linear flow rule l|3.10|) results in evolu- 
tion equations for F, G, and H: 

d t F ij = n 2k F kj + n jk F kl - u k Fj (5.3a) 

d t G l = fl lk G k - u k G l k + F l ^VL k l + (f V?*) (Q kj - u k ) (5.3b) 



d t H= 2n lk G k -u k H 



k+2{?tefy(Sl ki -u k j). (5.3c) 



Flow Rules. For compressible flows, two flow rules will be considered. We define 
them first, and then go on to consider their relative merits and demerits: 
I. Advection: ST 3 = 

II. Rotation: fi** = \ (u^ - u°^j 

Advection. For our anisotropic model, we shall advect £' and treat the quantities 
F, G, and H as parameters in the final system, each of which will have its own 
evolution equation. Substituting fi^ = u l - into the system Q5.3fl gives 

d t F = -£ U F (5.4a) 
d t G = -£ U G + F ■ grad(divu) (5.4b) 
d t H = 2 grad(div u) ■ G - u ■ grad H, (5.4c) 

One advantage of the advection flow rule is that it automatically closes the system 
(15.31) . For a general choice of f2, the system involves ^^C^k) anc ^ ^£'ffc£'^V which 
cannot be expressed solely in terms of F, G, and H . 

Rotation. For our isotropic model, we want to know whether the evolution equa- 
tion (|5.3a() for F preserves the isotropy relationship F — Identity. Suppose F u = 8^ 
at t = 0. Then substituting into Ij5.3a|l reveals that 

d t \t= F ij = Q ij + OP*. (5.5) 

If £1 is antisymmetric, we have dt\t=oF — 0, and F(x,t) = Identity solves l|5.3a|l for 
all t. We wish to know whether this solution is unique. This is guaranteed by a 
straightforward generaliza tion of the r esults concerning linear hyperbolic systems of 
first-order equations from lEvansI 1998j . assuming sufficient smoothness of u. 

We conclude that antisymmetry of fl is sufficient to guarantee that the initial 
data F — Identity is in fact preserved for all t. Then an immediate choice of a tensor 
f2 that is antisymmetric is given by the rotation flow rule pill . This form has a 
very attractive physical interpretation. Putting the linear flow rule equation l|3.10|l 
together with l|TTj) gives us 

where lu = curlu is the vorticity vector. The last equation can be interpreted in the 
sense that fluctuations are rigidly transported by the mean flow, with a local angular 
velocity given by the vorticity vector. 

Finally, the rotation flow rule lHT)l docs not by itself close the system (|5.3|) . When 
using this flow rule, we shall assume that G = and H — (3 2 . 

6. Equations for Averaged Dynamics. Here we shall write down two systems 
of coupled PDEs which describe the evolution of the average velocity and density in 
a compressible flow. Each PDE is derived from an associated averaged Lagrangian. 
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Compressible, fully Anisotropic, Inhomogeneous Fluids. By substituting (|TJ into 
the Lagrangian (|4.5|) , we obtain closure: the Lagrangian no longer depends explicitly 
on but instead on the tensors F, G, and H, for which a self-contained system 
of evolution equations l|5.4|l has already been derived — see 33 for details. Applying 
Hamilton's principle directly to l|5.2[) yields an evolution equation for u, the average 
fluid velocity. We write this equation using the operator A, which is defined as 



1 



(Av) 1 = - ( P v^ k ) 



(6.1) 



We also write w — pw' '(p) where ' means d/ dp as usual. The anisotropic compressible 
LAE-a equations are: 



-pw., 



a 

Y 



+ F$p ti w d + (*>) .. w + 2G%pw ti + 2G* {pw tn ). + {Hp 2 w') n \ (6.2a) 



d t p = - div(pu) 

d t F = -WF -u + F-Vu + Vu T ■ F 
d t G = -u-VG + G -Vu + F ■ grad(div <, 
d t H = 2 grad(div u) ■ G — u ■ grad H. 



(6.2b) 
(6.2c) 
(6. 2d) 
(6.2e) 



Well-posedness. We now sketch a rough well-poscdness argument for the system 
(16.21) . Assume that the tensor F is positive-definite. By this it is meant, since F 
is a (2,0) tensor, that for any one-form 8, the contraction F : (9 eg) 9) is positive 
everywhere. Given the p- weighted inner product (f,g)= J f g p, we have (/, —Af) = 
~ / / p {pf,jF jk ) k P ~ / f.]F jk f,kP > 0. Since —A is a positive definite linear 
operator, (1 — a 2 A) has trivial kernel and we expect that (|6.2|l is well-posed. 

It would be of analytical interest to see to what extent the "geode s ic par t" of 
these equations define a smooth spray in the sense of lEbin and Marsdenl [1970j, and 
which holds for the EPDiff equations (that is, the n-dimensional CH equations), as 
explained in lHolm and Marsdenl 2004l| . 

Compressible, isotropic, inhomogeneous. For this case we use flow rule (fTTf) , which 
can be written in vector notation as 



n 



Recall that this flow rule is compatible with an isotropic choice of the covariance 
tensor, i.e. F %3 — S 1 ^ . We further assume that G = and H = (3 2 for some constant 
(3. Using flow rule (fTTf) along with these extra assumptions in the general Lagrangian 
expression (15.21) gives us a Lagrangian in only two variables: 



l(u,p) 



A I 



-p\\u\\ 2 - pW(p) + a 2 



- / o(||VM|| 2 -tr(Vu- Vtt)) 



- \w'{p){\\Vp\\ 2 + p 2 p 2 )- l -w(p)^p 



d N x, (6.3) 
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where w(p) = W(p) + pW'(p) is the enthalpy introduced in l|4.6|) . Regarding this as 
a Lagrangian in u and u = p d N x, one uses the semidirect product Euler-Poincare 
equations fsee Ho lm. Marsden. and Ratiul 1998b]') to derive the system 

dt(pv) + (u ■ V) (pv) + a 2 div (pf7 • Vu) + pvdivu = — Vp (6.4a) 

d t p + div(pw) = 0. (6.4b) 

with the modified momentum pv and modified pressure p given by 

pv = pu + a 2 div (pSl) , (6-5) 

Vp = Vp + a 2 (3 2 pX7 ^ow' + ipV'^) . (6.6) 

Here are explicit coordinate expressions for two slightly complicated objects: 

pv 1 = pu l + ia 2 (p (u 3 z - u 1 ^ 
div (pfi • V«) = (/^Xj) k . 

The following convention for divergences of tensors has been used: given a 2-tensor 
, we set 

(dWAy =A% 

That is, the contraction implicit in the divergence operation always takes place on 
the first index. 
Observations. 

• In the case of homogeneous incompressible flow, where p is constant and 
div it = 0, the definition of pv in (|6.5(l reduces to 

v=(l- \a 2 ^j u, 

which after rescaling a to get rid of the factor of 1/2 is precisely the v one 
finds in treatments of the incompressible LAE-a and LANS-a equations. 

• The above does not work in one spatial dimension. The problem is that 
here fl reduces to (u x — u x )/2 — 0, which clearly does not describe transport 
at all. For a 1-D isotropic model one may very well want to forget about 
antisymmetry of il and instead use something such as the advection flow 
rule. One may, quite reasonably, conclude that the only meaning of isotropy 
in 1-D should be reflection symmetry. 

7. Future Directions. 

The Initialization Problem. Perhaps the largest unsolved problem for the La- 
grangian averaged equations is the initialization problem. A concise statement of the 
problem reads: 

Given initial data uq(x) for the Euler equation, how does one obtain 

initial data Uq(x) for the LAE-a equation? 
Let us look at this problem in slightly more detail. Let u denote the solution of the 
incompressible Euler equations for initial data uq, i.e. u(x,0) = uq(x). Similarly, 
let U denote the solution of the incompressible, isotropic LAE-a equations (|3.21|) for 
initial data Uq. 
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Now U should be, in some sense, the mean flow of the fluid. This means that 
Uq should be the mean flow of the fluid at time t = 0, implying that Uq should be, 
in some sense, an "averaged" or "filtered" version of uq. The question is: how does 
one derive Uq from uq! Another way of phrasing this question is: how do we describe 
(approximately) the initial state of the fluid (given exactly, for our purposes, by the 
field uq) using only the mean flow variable Uq? 

Numerous methods have been used to initialize the LAE-a equations for use in 
numerical simulations, but none of these methods has any theoretical foundation. 
There is also no theory regarding how one should filter a full Eulcr flow u, or even a 
family of flows u f -, in order to obtain a mean flow that could be compared with the 
full LAE-a trajectory U. In this respect, equation Q2.4JI . which states that 

(u £ o^(x,t))=u(x,t), 

is not helpful: we have no way to compute the fluctuation diffeomorphism group £ e . 
Therefore we have no way to compute the left-hand side (u £ o £ e ). 

The difficulty can be summarized in the following commutative diagram. Here S 
is the standard Euler action and S a is the Lagrangian-averaged action. 

Lagrangian average 
_ ga 

derive PDE, 
solve numerically 

the missing link 

Solid arrows represent steps that we know how to carry out. The dashed arrow 
represents the one step that we do not know how to carry out. Our strategy for this 
problem will be to develop methods by which we can test different filters for obtaining 
Uq from Mo in practice. 

Treatment of Densities. Another area for further investigation involves our treat- 
ment of the density tube p, e . There are two questions to ground us: 

1. We have tacitly assumed that at t = 0, and for all e, all x, 

H(x,t) = H e {x,t). 

An argument similar to the one made above in our discussion of the initializa- 
tion problem can be made here. Namely, p,(x, 0) represents the mean density 
at time t — 0. Meanwhile, p e (x,0) represents the true density of the fluid 
at time t = 0. These two quantities need not be equal. This prompts the 
question: how would we carry out the procedure from Sections |21 and ^ with 
tubes in which each trajectory does not have the same initial density p{x, 0)? 

2. As our derivation of the averaged compressible equations stand, we have 
derived the fact that the "mean" density \i was advected by the mean flow U: 
dtfi = —£u(i. Substituting /j, = pd N x and using the definition of divergence 
yields the standard continuity equation 

dtp + div(pU) = 0. 

In both RANS and LES treatments of averaged/filtered flow, the mean flow U 
satisfies a modified continuity equation rather than the standard one. There- 
fore: why does the Lagrangian averaged mean density p, satisfy the usual 
continuity equation? 



S 
u 
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The two questions regarding densities are in fact related. To see this, let us suppose 
that given the initial density po associated with the center line of our tube rj, we have 
a method for constructing a family of initial densities Pq for each of the other curves 
in the tube rf . Now defining 4 

M £ (t) = (% e )*Mo and m = ( M e (t)), 
we will find that p,(t) = p{t) d N x satisfies a modified continuity equation 

d t p(t) + div(pu) + div ^p e (eu 1 + ~eV'J J = 0. 

To close this equation, we must either carry out the average directly, or we must 
expand p e about a suitable trajectory and make modeling assumptions. 

Filtered Lagrangians. We have seen that the current averaging procedure leads to 
complicated averaged equations. Furthermore, there is no clear way to evaluate nu- 
merically the flow rules we have proposed on physical grounds. One of our immediate 
goals is to investigate a filtering approach, still at the level of the Lagrangian, which 
will lead to simpler averaged models that can be tested numerically. The filtering 
approach we have in mind begins with a decomposition of the velocity field 

u — u + u' and p = p + p' (7.1) 

into mean and fluctuating components. This would replace the Taylor expansion 
(|4.2[1 of u e and p e that we carried out in the present work, and would therefore lead 
to Lagrangians and equations with much less algebraic complexity. As opposed to the 
axiomatic averaging operation (•), the filter shall be specified concretely. We expect 
this to help greatly with the initialization and density problems discussed above; 
furthermore, the filtering approach leads naturally to questions about the relationship 
between LES and LAE-a models. 

Simpler Models. As we previously noted, the flow rule approach developed in 
this paper does not yield a one-dimensional compressible averaged model. We are 
currently investigating such a model, derived from the filtered Lagrangian 

l(p,u) = J (^uv~W(p)^pd N x, (7.2) 

where v = (l — a 2 d xx )u. To derive this Lagrangian, we filter only the velocity, leaving 
density and potential energy alone. This is the compress ible analogue of the filtered 
Lagrangian used in deriving the Camassa-Holm equation ICamassa and Hold |1993| . 
The analysis and numerical simulation of the new equations presented in Section 
of this work will be difficult. Much easier is the analysis of the PDE associated with 
(|7.2(l . In particular, we expect that numerical studies of this one-dimensional model 
will yield insight into the dynamics of the higher-dimensional equations. 

Entropy. In the derivation of our compressible averaged models, we have made 
the barotropic assumption W = W(p). We expect the resulting barotropic model to 
be useful in computing mean flow quantities in regimes where we are not concerned 
with strong physical shocks, for example in climate models. The next major step 
forward will be to remove the barotropic assumption, and derive a model that is valid 
in regimes where we are concerned with shocks. 



"Note that (^{x,t)) ^ fJ.(x,t). 
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To this end, we have derived an averaged model for the general case, where 
the potential energy has the form W(p,S), where S is the entropy. This model, 
which consists of a system of equations for p, u, and S 1 , also involves the pressure 
p. Therefore, in order to close the system, we require an equation of state relating 
p to p and S. The open question now is as follows: given an equation of state for 
the compressible Euler system, what is the equation of state relating the averaged 
variables to one another? In other words, how does Lagrangian averaging interact 
with the thermodynamics of the system? We hope that analyzing a finite-dimensional 
case of this interaction will shed light on this issue. 

Connections with Kevrekidis ' Coarse/Fine Methods. Given a description of any 
mechanical system, not necessarily involving fluids, in the form of a Lagrangian £, we 
can carry out the procedure described in |J5]to find an averaged Lagrangian (I). From 
this we can derive equations of motion for the average dynamics of the original system. 
Changing our language slightly, we say that we have a general method for extracting 
the "coarse" dynamics of a mechanical system whose full description involves motions 
on both fine and coarse scales. 

Another method fo r computing the coarse-s cale dynamics of a mechanical system 
has been put forth by in lKevrekidis et. al.l|2003| . Kevrekidis' method does not involve 
trying to write down equations of motion which govern the coarse dynamics. Instead, 
he offers an algorithmic approach, the crux of which is as follows. The coarse dynamics 
of a system are found by lifting the initial (t — to) state to an ensemble of initial states, 
integrating each using the full equations until some small final time t = e has been 
reached, and projecting the resulting t = e states onto a single state. This t = e state 
is then extrapolated to a state at some desired t = tf > 0. By iterating this process 
and tuning the lifting, projection, and extrapolation operations, this method can be 
used to recover the coarse dynamics of the system. 

Now the question that begs to be asked is as follows: for the case of fluid dynamics, 
how different are the coarse dynamics provided by the LANS-a equation from the 
coarse dynamics one would obtain by following Kevrekidis? The difficulty in answering 
this question lies in implementing a full fine-scale integrator for fluids that one could 
successfully embed inside Kevrekidis' coarse-scale algorithm. We look forward to 
tackling this task soon. 
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9. Appendix: Fluctuation Calculus Details. Before proceeding with any 
derivations, we state the Lie derivative theorem for time-dependent vector fields: if 
the vector field X\ has flow F\, then 



Our task now is to derive equations (|2.10(1 . Starting with 1)2. 6|l , let us move rf to 
the right-hand side of the equation: 




(9.1) 



(9.2) 
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The strategy is to differentiate with respect to e and use the Lie derivative theorem 
(|9.1|) . The intrinsic definition of divergence 

£ c (y) = (div„C)i> ( 9 - 3 ) 

and the canonical volume form d X — dx 1 A - ■ -Adx N will both be used in what follows. 
Note that div£ with no subscript on the div means £((d N x). Before applying the Lie 
derivative theorem, note that the vector field 

W = ^ (,e)-i (9.4) 



has flow rf . A simple computation yields 

W e = £" - V£' • (9-5) 



0_ 

7h 



e=0 



Then we start with p': 



d d 

— fj,Q = = — (?7 C )* p e by differentiating (|9.2j) 

= (vr(j£r + £w<^ by ED 

= ?y* (/t + £^ijS) at e = 
=> fjf = —£^p 

p 1 d N x = - (£ e p) d N x - p (£ e d N x) by (£3 

// d^x = - (Vp ■ £' + p div £') d^a; by 

=>y = -divoo 



Next we compute p' 

^A* = ° = £5 Of) A*' 



0?T ( ^tA* £ + £w^- + ^- {£w^) + £w£wfJ> ( 



v <9e 2 de de 

=> = rj* (//' + 2£%p! + .££/'_v£'.f'|U + £^£^>p) 
=> //' = —£^/fi + 2£^/£^p — £^i£^>p + £y£r.£i/j, 
p" d N x = - (div (pf)) d N x + div (div (p?) C) d N x + div (pV£' • £') ^ 



-(div(pn) + ((pr).e ,J j • (/<rr) 
-(div(/«")) + (^V) 



= - (div (p?')) + div div (p? ® O 
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